%% 对混频信号和滤波信号做时域分析
clc                 % 清屏
clear all;          % 清除变量
close all;          % 关闭所有窗口
% 加载文件
load before_filter_data.txt
load after_filter_data.txt
% 给文件命名
before=before_filter_data;
after=after_filter_data;
%% 绘图
figure(1);
plot(before);                       % 绘制混频信号
axis([0 200 -20000 20000])          % 坐标轴范围设置
title('滤波前混频信号时域波形');      % 标题

figure(2);
plot(after);                        % 绘制滤波后信号
axis([0 200 -2*10e8 2*10e8])        % 坐标轴范围设置
title('滤波后信号时域波形');          % 标题

%% FFT分析
fs=50e6;                % 采样频率
len=12255;              % 滤波信号长度
len2=12255;             % 混频信号长度
ff=fs/len;              % 滤波信号频谱分辨率
ff2=fs/len2;            % 混频信号频谱分辨率
f=-fs/2:ff:fs/2-ff;     % 滤波信号频谱范围
f2=-fs/2:ff2:fs/2-ff2;  % 混频信号频谱范围

%% FFT
bfft=fft(before);       % 混频信号做FFT
afft=fft(after);        % 滤波信号做FFT
figure(3);
plot(f2,abs(fftshift(bfft)));   % 绘制混频信号频谱
axis([-8e6,8e6,-inf,inf]);      % 坐标轴范围设置
title('滤波前混频信号频谱');      % 标题
figure(4);
plot(f,abs(fftshift(afft)));    % 绘制滤波信号频谱
axis([-2e6,2e6,-inf,inf]);      % 坐标轴范围设置
title('滤波后单频信号频谱');      % 标题